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I. INTRODUCTION 

Understanding the internal structure of nucleons has been a central research topic in nuclear and particle physics 
for many decades. As detailed experimental data continue to emerge, improved theoretical understanding of the 
hadronic spectrum will be needed to learn more about the complex, confining dynamics of quantum chromodynamics 
(QCD). Lattice calculations offer a means of linking experimental data to the Lagrangian of QCD, allowing access to 
the internal structure of any resonance. 

At Jefferson Laboratory ( JLab) , an ambitious program of research into a range of hadronic excitations is underway. 
To date, the Hall B experiment has collected a large amount of data regarding the spectrum of excitations of the 
nucleons. The Excited Baryon Analysis Center (EBAC) [3, Q aims to review all observed nucleon excitations system- 
atically and to extract reliable parameters describing transitions between resonances and the ground-state nucleons. 
The 12-GeV upgrade of JLab's CEBAF accelerator will make possible the GlueX experiment, which will produce an 
unprecedented meson data-set through photoproduction. A particular focus will be the spectrum of hadrons with 
exotic quantum numbers which can arise when the gluonic field within a meson carries non- vacuum quantum numbers. 
Such "hybrid" mesons offer a window into the confinement mechanism and will be studied theoretically in some detail 
using lattice methods. Lattice spectroscopy can determine the properties of isoscalar mesons as well, including any 
possible candidate glueballs. 

Accurate resolution of excited states using lattice QCD has proven difficult. In Euclidean space, excited-state 
correlation functions decay faster than the ground-state and at large times are swamped by the larger signals of 
lower states. To improve the chances of extracting excited states, better temporal resolution of correlation functions 
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is extremely helpful. An anisotropic lattice, where the temporal domain is discrctized with a finer grid spacing 
than its spatial counterpart, is one means of providing this resolution while avoiding the increase in computational 
cost that would come from reducing the spacing in all directions. This finer resolution must be combined with 
application of variational techniques to construct operators that overlap predominantly with excited states. In a 
series of papers 0, 0, HI, H, Q, techniques to construct operators in irreducible representations of the cubic group to 
extract radially and orbitally excited states were presented. Application of these techniques to quenched anisotropic 
lattices shows clear signals as high as the eighth excited state. Also possible on the anisotropic lattice are studies of 
radiative transitions in meson systems^ and nucleon-Pn transition form factors^. 

Making the lattice discretization anisotropic comes with a price, however. Since hypercubic symmetry is broken 
down to just the cubic group, relevant (dimension-four) operators can mix in the lattice action. To ensure the 
continuum limit of the lattice theory has full Lorentz invariance, a nonperturbative determination of the lattice action 
parameters that enforce the symmetry at finite lattice spacing in some low-energy observables has been performed [IcI]. 
In work to be reported elsewhere, a perturbative determination of these action parameters is also being carried out 
by this collaboration . 

In this study, we perform three-flavor dynamical calculations with two degenerate light quarks and a strange quark. 
In a previous study, we tuned a three-flavor lattice action to ensure Lorentz symmetry is restored in appropriately 
chosen low-energy observables. We showed empirically that restoring the symmetry at quark masses below 175 MeV 
requires only small changes to the action parameters, and no further determinations of these parameters is needed 
within the scope of this study. Our fermion action is a Sheikholcslami-Wohlert discretization, generalized to the 
anisotropic lattice[l2|]. The fermion fields interact with the gluons via 3-dimensionally stout-smeared [HJ links, and 
the gluon action is Symanzik improved at the tree level of perturbation theory. To assess the cost of these dynamical 
calculations, we study the efficiency of the hybrid Monte Carlo (HMC) algorithm in large-scale production simulations. 

Using current algorithms and computing resources, it remains impractical to run calculations at the physical value of 
the light-quark mass. An extrapolation of the light-quark dependence of simulation data is needed. While simulations 
straddling the correct strange-quark mass have been performed, determining the appropriate choice of mass in the 
Lagrangian is also problematic; a priori, this value is not known and changing the bare strange-quark mass affects 
all lattice observables in a delicate way. This work proposes a simple means of setting the lattice strange-quark 
mass by examining dimensionless ratios with mild behavior in the light-quark chiral limit. A new set of coordinates, 
parameterizing the space of theories with different light and strange-quark masses is introduced to help this process. 
We use the ratios Iq — 9rn^/4rriQ and sq = 9(2m K — mfy/Am^, inspired by expanding the pseudoscalar meson masses 
to leading order in chiral perturbation theory. 

With this framework in place, the spectrum of some ground-state mesons and baryons is determined and extrapo- 
lated to the physical quark masses using leading-order chiral perturbation theory. The Sommer scale 1J| is determined 
(in units of the Omega-baryon mass) on a subset of our ensembles and extrapolated to the physical quark masses. 
Using this method, we intend to continue our exploration of excited-state hadrons, including isoscalar and hybrid 
mesons. A clear means of handling unstable states is needed, and our suite of measurement technology is currently 
under further development. A study of these techniques on Nf = 2 dynamical lattices gives us confidence that more 
precise understanding of these states will be forthcoming. 

The structure of this paper is as follows: In Sec. [Ill wc wm discuss the actions and algorithms used in this work, and 
the performance of the method used to generate Monte Carlo ensembles is examined. The details of the measurements 
we performed on these ensembles is presented in Sec. IIIII Sec. II VI presents the method we propose to set the strange- 
quark mass, including the dimensionless coordinates used for extrapolating quantities measured at unphysical quark 
masses to the physical theory. Our determination of a selection of states in the hadron spectrum and the Sommer 
scale is given in Sec.[V] Some conclusions and future outlook are presented in Sec. I VII 



II. SIMULATION DETAILS 



In this section details of the lattice action and the performance of the hybrid Monte Carlo algorithm are presented. 
Monte Carlo simulations were performed on lattices with grid spacings of a s and at in the spatial and temporal 
directions respectively and with physical volumes L 3 x L t where L s = N s a s and L t = N t a t . Lattices with extents 
xN t = 12 3 x 96, 16 3 x 96, 16 3 x 128 and 24 3 x 128 were employed. 



Action 



The gauge and fermion actions used in this work are described in great detail in our previous work [1(1. For 
completeness in this paper, we briefly review the essential definitions. For more detailed definitions, see Ref. jlOj ] - 
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For the gauge sector, we use a Symanzik-improved action with tree-level tadpole-improved coefficients: 
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where fivy = ReTr(l — W) and W = V, the plaquette, or 7^,,, the 2x1 rectangular Wilson loop (length two in 
the fj, direction and one in the v direction) with {s,s'} G {x,y, z}. The parameter j g is the bare gauge anisotropy, 
N c = 3 indicates the number of colors, j3 is related to the coupling g 2 through j3 = 2N c /g 2 , and u s and ut are the 
spatial and temporal tadpole factors. This action has leading discretization error at 0{a^, a 2 , g 2 a 2 s ) and possesses a 
positive-definite transfer matrix, since there is no length-two rectangle in time. 
In the fcrmion sector, we adopt the anisotropic clover fermion action (l2|: 
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1 1 

7/ «2 f£, 



(2) 



where jf is the bare fermion anisotropy and £r = a s /at is the renormalized anisotropy. J s ,t-> a st and a ss : (with 
a iiv = ^[7/ii7i/]) are Dirac matrices. Hats denote dimensionless variables which connect to dimensionful quantities 



as: quark field tp 
"Wilson operator" 
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ip, bare quark mass mo : 
V„, - |7^A M (with V M 



moat, gauge field strength F M „ = a ll a u F llu = jluv{V^ u {x)) and 
: a^V^, A M = a^A^). The gauge links in the fermion action are 



3-dimcnsionally stout-link smeared gauge fields with smearing weight p = 0.14 and n p = 2 iterations. u s and u t are 
the spatial and temporal tadpole factors from smeared fields, respectively. 
In our previous work[l(|, we found at /3 = 1.5, that the tadpole factors are 



u s = 0.7336, u t = 1, u s = 0.9267, u< = 1. 
Tuning the anisotropy for all quark masses (even below the chiral limit) gives the desired 7* * 

7 a *=4.3, 7; = 3.4. 



(3) 
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B. Algorithm 

We use the rational HMC (RHMC) algorithm for gauge generation [l^. The theoretical aspects of our procedure 
were discussed in detail in Ref. [l(j. Here we discuss only the aspects that are specific to the calculations presented 
in this work. 

We use rational approximations for both the light-quark fields and for the strange quarks — one field for each light- 
quark flavor and another one for the strange. We employ even-odd preconditioning for the Wilson clover operator, 
obtaining the Hamiltonian 

# = ^Tr7rt7r-2^Tr log A ee (mi ) - ^ Tr log A ee (m s ) + S P (mi) + S F {mi) + S F {m s ) - S S G - (5) 

X^L X X 

where tt are the momenta conjugate to the gauge fields; terms involving A ee contribute effects due to the parts of the 
preconditioned clover determinant coming from the submatrix connecting even sites; Sq and S G are the parts of the 
gauge action involving loops in the spatial directions only and with loops including time direction respectively; and 
Sf{i^i) and £p(m. s ) are pseudofermion terms for the rational approximations to the fermion action corresponding to 
the light- and strange-quark fields respectively, which we discuss below. 

The pseudofermion terms Si?(m) employ a rational approximation to the fermion determinant of the even-odd 
preconditioned clover operator coming from the submatrix connecting the odd sites for a quark with mass m. This 
submatrix is 

M(m; U) = A 00 (m; U) - D oe {U)A^{m- U)D eo (U), (6) 

where A oa is the clover operator on the odd sites, A~^ is the inverse clover operator on the even sites and D oe (D eo ) 
is the Wilson hopping term connecting odd sites with even (even sites with odd). In all the expressions, U denotes 
stout-smeared gauge fields U, which were smeared as described in Subsection III Al 
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V 


light quark 
atmi bounds no. poles max. error 


strange quark 
atm 3 bounds no. poles max. error 


24 3 x 128 


-0.0840 (5 x 10~ 6 , 10) 16 1.8 x 10~ 8 


-0.0743 (10^ 4 , 10) 12 8.8 x 10~ 8 



TABLE I: Details ol the approximations used for the pseudofermionic action R * . We show the bounds, the number ol poles 
and the maximum error for the approximation to the light and strange pseudofermion terms. 



V 


light quark 
atmi bounds no. poles max. error 


strange quark 
atm 3 bounds no. poles max. error 


24 3 x 128 


-0.0840 (5 x 10~ 6 , 10) 12 2.5 x 10~ 6 


-0.0743 (10~ 4 , 10) 10 2.0 x 10~ 6 



TABLE II: Details of the approximations used for the force R 2 . We show the bounds, the number of poles and the maximum 
error for the approximation to the light and strange pseudofermion terms. 

To construct our pseudofermion actions, we use the rational approximation Rt(M'M) in partial-fraction form: 

R% (M f M) = a Y^Pi {M^M + q^ 1 w (Aft M ) f , (7) 



where wc drop the quark-mass dependence of M for clarity. The coefficients a, pi and qi define the approximation and 
are determined via the Rcmez algorithm (l6l. \v\ applied over the spectral bounds of the operator M^M. In particular, 
we needed to compute approximations with (a, b) = (—1,4) for evaluating the actions (see below), (a, b) = (1,4) for 
pseudofermion refreshment and (a, b) = (—1, 2) for our molecular dynamics (MD). Our approximation bounds for the 
action are shown in Table |TJ We solve the linear system resulting from applying Rt to pseudofermion fields using the 
multi-shift conjugate gradient algorithm [3]. We use a stopping relative residuum r < 1CP 8 in our energy calculations 
where the residuum for pole i is 



\<t> - (M+M + qi) i>i 



(8) 



where <p is the pseudofermion field and ipi is the solution corresponding to the i pole. However, since the multi-shift 
algorithm cannot be restarted, our stopping was based on estimates of accumulated with the short-term recurrence 
in the solver algorithm, which may be slightly different from the true residual as defined in Eq. [8] due to solver 
stagnation and rounding effects. To minimize rounding effects we accumulated sums and inner products using double 
precision. 

Our pseudofermion action terms are 



S F = X^X, X = R-i (APM) <f>= a^piijji, 



(9) 



individually for each flavor. We do not need to employ multiple pseudofermion fields per flavor in this study. 

During our simulation, we adjusted our approximation range by measuring eigenvalue bounds every five trajectories 
during the process of thermalization. Thereafter we continued to measure the bounds to ensure we do not suffer from 
boundary violations. 

Our molecular dynamics process employs a rational force 



F = -a^pitpl ( 



dM 1 



■M 



(10) 



where i runs over the number of poles in the approximation _R~2. Time derivatives are evaluated over the stouted 
gauge field U, and only the final sum is recursed down to compute the force for the thin links U. 

We employ a multiple-timescale integration scheme for the molecular dynamics evolution [l9j by nesting a second- 
order Omelyan[2(| [2l| integration step at each timescale. Our largest forces come from the temporal directions: 
the gauge force from Sq and the temporal forces generated by the pseudofermions. To mitigate the numerical effort 
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V 


(a t mi,a t m s ) 


poles for atmi 


residua 


poles for atm 3 


residua 


24 3 x 128 


(-0.0840, -0.0743) 


12 


10- 4 ,10 _4 ,5 x 10" 5 , 
5 x 10~ 5 ,5 x 10" 5 , 10" 5 , 
10~ 5 ,5 x 10" 6 ,5 x 10" 6 , 
5 x 10" 6 ,3 x 10" 6 , 10" 6 


10 


10" 4 ,10 _4 ,5 x 10 -5 , 
10" 5 ,10- 5 ,10- 5 , 
5 x 10 _6 ,5 x 10 _6 ,3 x 10" 6 , 
10" 6 



TABLE III: Requested residua for the poles in the MD force approximation from the smallest shifts (leftmost) to larger shifts 
(rightmost). 



(atmi y atm a ) dt\ dt 2 /dt\ £md Accept Rate 



V 



24 3 x 128 (-0.0840, -0.0743) 



3.5 



0.71 



TABLE IV: The two timescales used in the molecular dynamics integration. The spatial timestep for the coarse scale is dtl, 
and for the finer scale it is dt 2 s , which we display as a fraction of dtl here. We also show our MD timestep anisotropy. On each 
scale, dtl = dt\ /£md- Finally, we show the average acceptance rate for the molecular dynamics with these step sizes. 



needed 22|, we place the S G term in the action on a finer timescale than the other terms, and to deal with the temporal 
forces from the pscudofcrmion terms, employ an anisotropic timestep with temporal timestep dt of length 



(11) 



where £md = 3.5. Apart from the above, we find the forces from Sq and the spatial forces from the Sf terms to be 
within a factor of 2 of each other, so we place them on the same timescale. The forces from the Tr log A ee terms are 
very small in comparison but have small numerical cost, so we place them on the same timescale as S G and Sf- 

We note that the Hamiltonian for the MD does not need to be known as accurately as the one for the energy 
calculations; all that is required is for the MD to be reversible, area-preserving, and (as a practical matter) for the 
acceptance rate to be reasonable. To save on numerical effort we solved our systems of linear equations only to a 
residuum tmd of at most tmd < 10~ 6 . Correspondingly, we never required the rational approximation to R~? to 
have a maximum error better than 10 -6 , resulting in a smaller number of poles in the approximation than we need for 
the energy calculations. Further, to make the MD even less numerically intensive, we follow Refs. [H, [24[ by relaxing 
the requirements on the residua for individual poles in R~? . We use a range of tmd < 10 -4 for the smallest shifts 
and tmd < 1CP 6 for the larger shifts. We tune our molecular dynamics to attain an overall acceptance rate close 
to 70%. We show in Table ITT1 the bounds of the MD rational approximation used for the run with V = 24 3 x 128, 
(a t mi,a t m s ) = (—0.0840,-0.0743). We show the residua requested in the MD evolution in Table lllll and the timesteps 
and the resulting acceptance rate in Table ITVl 



C. Thermalization and Autocorrelation 



During the first segment of each gauge ensemble generation, some special conditions apply. We do not apply the 
acceptance test during the first O(10) trajectories in each series, which allows a fast initial approach to the vicinity 
of the equilibrium. Such a scheme is particularly important in the case of simulations starting from totally ordered 
or disordered configurations. Wherever possible, however, we begin the algorithm with an equilibrated configuration 
from a simulation at nearby parameters. Also, during this phase (as mentioned above), the minimum and maximum 
eigenvalue bounds are updated every 5 trajectories. 

Figure [1] shows its plaquette history for 24 3 x 128 volume and atirii = —0.0840. Both plaquette histories (one 
excluding temporal links and the other including only plaquettes with temporal links) show that equilibrium is 
reached long before 1000 RHMC trajectories. Therefore, to allow for thermalization of our gauge ensembles during 
the RHMC, we discard the initial 1000 trajectories from each set. 

Figure [5] shows a histogram of the lowest eigenvalues of the Dirac operator M^M for the light and strange quarks 
from the ensemble with 24 3 x 128 volume and atvni = —0.0840. The lowest eigenvalues remain safely above the 
minimum eigenvalue bounds in which our rational approximation is valid. In addition, they show a clear gap away 
from zero, where the stability of the algorithm might be compromised. 
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FIG. 1: Temporal (left column) and spatial (right column) plaquette history from the ensemble with 24 J x 128 volume and 
atrni — —0.0840. The x-axis is in units of trajectories. 
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FIG. 2: Histogram of the lowest eigenvalues of the Dirac operator for the light and strange quarks from the ensemble with 



24 x 128 volume and a t mi 



-0.0840. 



The autocorrelation function is defined as 

P {t) = {(0(t')-{0))(0(t' + t)-(0))}, (12) 

where (...) means taking an average over the samples, t is the trajectory difference in the autocorrelation (from 1 
to TV total trajectories), and different t' (also indexing trajectory number) are averaged. To calculate the integrated 
autocorrelation length T; nt with jackknife-cstimatcd errorbar, we first divide the configurations into blocks of size TV;,; 
we calculate Pj(t) for jackknife index j by ignoring contributions when either t or t' + 1 is located within the j th block 
and replacing (O) by (0)j, the mean value without the j th block. With a jackknife data set of length TV/TV;,, we 
calculate integrated autocorrelation length, 

r int (i max ) = 5 + p (f = o) ( 13 ) 

using standard jackknife procedure. The autocorrelations of the spatial plaquette from gauge ensemble atmi = 
—0.0808, 16 3 x 128 are shown in Figure [3l the integrated autocorrelation length for the stout-smeared plaquette is 
about 30 trajectories, which is around twice as large as the un-smeared ones. The integrated autocorrelation length 
for lowest light and strange eigenvalues are around 13 and 10 trajectories respectively; shown in Figure [H Figured] 
shows the case of pion and proton correlators at t = 30 on our largest spectrum measurement (518 configurations) 
ensemble, a t mi — —0.0808, 16 3 x 128. The integrated autocorrelation length is about 30 trajectories. 




FIG. 4: Autocorrelation p(t) and integrated autocorrelation length Tint (in trajectories) for the up/down (above) and strange 
(below) quark eigenvalues from the ensemble with f6 3 x f28 volume and a t mi — —0.0808. 




III. MEASUREMENT TECHNIQUES 



In this section, the methods used to determine relevant spectroscopy data on the Monte Carlo ensembles are 
described. We have used well-established lattice spectroscopy technology throughout this calculation. 

To better access ground-state correlation functions, we use the variational methodj25l.l26l|. Consider the generalized 
eigenvalue problem 

C(t)v = \(t,t )C(t )v, (14) 

where to is chosen as the earliest time at which our model (given below) well describes the correlator C. (t) is a 
two-point correlation function, composed from the operators Oi and Oj. The correlation matrix can be approximated 
by a sum over the lowest N states: 

oo 

n=l 
N 

« Y, u ™*v Enit ~ to) > ( 16 ) 

n=l 

where E n is the energy of the n th state, and u m ■ z n = S mn . We extract the energies from the eigenvalues 

A„(Mo) = e-^'-'o), (17) 

which are obtained by solving 

C(t )- 1/2 C(t)C(t Q )-^ 2 v n = A„(t, to)v n . (18) 



A. Hadron Correlation Functions 



We perform measurements starting from trajectory 1000 on every 10 th trajectory, using the EigCG inverter (devel- 
oped by A. Stathopoulos et al. in Ref. [27]]) to calculate quark propagators (with CG residual set to 10~ 8 ). We use 
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FIG. 6: Pion and rho (above) and proton and Delta (below) fitted masses as functions of t„ 
t — 0. The bands indicated the final fitted masses summarized in Tables IVT1 and IVTT1 
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TABLE V: Meson interpolating operators. States are sorted into columns according to the degree of strangeness from (left 
two columns) to 1 (right column) and then according to total isospin. 



4 sources on each configuration, where a random source location is selected for the first source, and the remaining 
three are uniformly shifted by N XViZ /2 and N t /A\ this arrangement should reduce potential autocorrelations between 
configurations. We bin the data over spans of 5 measurements. 

In this work, we construct a 3 x 3 correlator matrix Cy by using 3 different Gaussian smearing widths (a G 
{3.0, 5.0, 6.5}) on the hadron operators. We extract the ground-state principal correlator and fit the ground-state mass 
using a cosh form. (We also try an exponential form on the principal correlator, and the fit results are consistent.) 
The i m i n dependences (with i max f=a 50) of the fitted masses are shown in Figure [5] for the pion, rho, nucleon and 
Delta. The fitted masses are very consistent between various choices of starting time in the fits. 

We use meson interpolating fields of the form qTq, which overlap with the physical states listed in Table fVl charge 

conjugation C applies only to particles with zero net flavor. The estimated i] mass is </ \m\ + ^m^g. The ground- 
state masses are summarized in Table IVT1 We have two volumes (12 3 and 16 3 ) of the lightest ensemble, m s = —0.0540, 
and two (16 3 and 24 3 ) on atm s = —0.0743: no major finite-volume effects are observed, except for the clq mass from 
the atm s = —0.0743 ensemble and baryon states from atm s = —0.0540. 

The octet baryons are calculated using the interpolating field ( qi 6*7475(72) q 1 (with qi = u/d or s quark); the A uses 
2(uCj5d)s+ (sCj5d)u+ {uC"f^s)d; and the decuplet uses 2(<7 2 C(1/2)(1 + 74)7-91 )qi + {qiC( 1/2) (1 +74)7- qi)q2 (with 
7_ = 7 X — j y ). The calculated octet and decuplet ground-state masses are summarized in Table Ivnl We observe 
a finite-volume discrepancy in the baryon sector on the lightest ensemble, atm s = —0.0540. When we extrapolate 
the hadron masses to the physical limit, we will exclude the small volume sets: 12 3 with atm s = —0.0540 and 16 3 
with atm s = —0.0743. Figures [7] and [8] show the squared-pion-mass dependence of these quantities. We note that for 
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U.ZOZZl_LO) U.ZOZZI±Ot U.ZOZZl±OI 


0.286(5) 


0.286(5) 0.286(5) 


0.415(20) 0.436(15) 0.459(18) 


81 9(1 9^ 


^0 


16 


128 


-0.0743 


-0.0743 


0.1483(2) 0.1483(2) 0.1483(2) 


0.2159(6) 


0.2159(6) 0.2159(6) 


0.287(6) 


0.317(5) 0.325(5) 


0.6867(17) 


79 


16 


128 


-0.0808 


-0.0743 


0.0996(6) 0.1149(6) 0.1196(5) 


0.173(2) 


0.1819(21) 0.1901(18) 0.222(11) 


0.252(6) 0.269(5) 


0.574(6) 


518 


16 


128 


-0.0830 


-0.0743 


0.0797(6) 0.1032(5) 0.1100(4) 


0.1623(16) 0.1733(10) 0.1845(11) 0.196(18) 


0.236(8) 0.263(8) 


0.491(6) 


266 


16 


128 


-0.0840 


-0.0743 


0.0691(6) 0.0970(5) 0.1047(5) 


0.154(3) 


0.1663(16) 0.1788(13) 0.159(15) 


0.222(7) 0.238(8) 


0.448(7) 


224 


24 


128 


-0.0840 


-0.0743 


0.0681(4) 0.0966(3) 0.1045(3) 


0.1529(10) 


0.1660(6) 0.1788(6) 


0.194(14) 


0.233(4) 0.242(6) 


0.446(3) 


287 



TABLE VI: Meson masses for Nf = 3 and Nf =2 + 1 (in temporal lattice units). 



N a 


N t 


a t mi 


a t m s 




a t m N 


a t m s 


a t m= 


a t m,A 


a t 7TlA 


a t m s * 


a t ms* 


a t rriQ 


^cfg 


12 


96 


-0.0540 


-0.0540 


0.2781(9) 


0.521(4) 


0.521(4) 


0.521(4) 


0.521(4) 


0.556(7) 


0.556(7) 


0.556(7) 


0.556(7) 


92 


12 


96 


-0.0699 


-0.0540 


0.1992(17) 


0.398(5) 


0.420(4) 


0.439(4) 


0.418(4) 


0.452(6) 


0.470(7) 


0.487(6) 


0.501(4) 


110 


12 


96 


-0.0794 


-0.0540 


0.1393(17) 


0.318(6) 


0.356(5) 


0.386(4) 


0.351(5) 


0.365(10) 


0.398(9) 


0.424(7) 


0.452(5) 


95 


12 


96 


-0.0826 


-0.0540 


0.1144(19) 


0.295(11) 


0.338(9) 


0.369(6) 


0.330(7) 


0.353(10) 0.382(11) 


0.414(8) 


0.447(5) 


84 


16 


96 


-0.0826 


-0.0540 


0.113(3) 


0.273(8) 


0.316(7) 


0.350(5) 


0.310(6) 


0.309(10) 


0.347(9) 


0.385(7) 


0.423(8) 


25 


12 


96 


-0.0618 


-0.0618 


0.2322(15) 


0.433(7) 


0.433(7) 


0.433(7) 


0.433(7) 


0.470(8) 


0.470(8) 


0.470(8) 


0.470(8) 


50 


16 


128 


-0.0743 


-0.0743 


0.1483(2) 


0.3165(18) 0.3165(18) 


0.3165(18) 0.3165(18) 


0.353(3) 


0.353(3) 


0.353(3) 


0.353(3) 


79 


16 


128 


-0.0808 


-0.0743 


0.0996(6) 


0.242(4) 


0.259(4) 


0.266(3) 


0.256(3) 


0.284(8) 


0.297(6) 


0.304(5) 


0.311(6) 


518 


16 


128 


-0.0830 


-0.0743 


0.0797(6) 


0.220(3) 


0.242(3) 


0.2510(19) 


0.236(2) 


0.270(7) 


0.283(5) 


0.292(4) 


0.304(3) 


266 


16 


128 


-0.0840 


-0.0743 


0.0691(6) 


0.207(4) 


0.229(3) 


0.239(2) 


0.218(3) 


0.262(8) 


0.275(7) 


0.282(5) 


0.294(4) 


224 


24 


128 


-0.0840 


-0.0743 


0.0681(4) 


0.2039(19) 0.2287(15) 


0.2395(12) 0.2209(15) 


0.256(3) 


0.271(3) 


0.282(2) 


0.2945(16) 


287 



TABLE VII: Baryon masses for Nf = 3 and Nf =2 + 1 (in temporal lattice units). 



non-strange hadrons, the sea strange-quark dependences are relatively mild. 

At low pion masses, not all the states we calculate on the lattice are safe from decays. To check which particles 
may decay, we compare the particle masses to the threshold two-particle energies in each channel. The vector mesons 
could decay to two pseudoscalars in a P-wave: p — > 7t(p) + tt(— p), K* — ► K(p) +7r(— p) and <fi — > K (p) + K(— p), where 
the minimum allowed momentum p on the lattice is In Figure [3 we plot the lowest two-particle energy threshold 
for the m s = —0.0540 data with our two lattice extents (dot-dashed line) and atm s — —0.0743 (dashed). All vector 
mesons in our calculation are well below threshold. The scalar mesons could decay to 7777 in an S'-wave, which puts 
the states slightly below the threshold. Similarly for the a\ and 61 mesons: a\ — > 7r(0) + p(Q) and 61 — > 7r(0) +u(Q). 
In this case, we approximate the u> by the p, since their masses are similar. The do, a\ and 61 (especially for the 
m s = —0.0743 ensemble) are slightly below the decay threshold. Fortunately, we have two volumes on the lightest 
mi set; if these states are not single-particle, the ratios of their overlap factors between the two volumes would be 
of order 2 or highcr[28]. We find the ratios (using point-point correlators) to be 0.87(12), 0.97(6) and 1.01(6), which 
indicates that our measurements are of single-particle states. The decuplet baryons are free from decays into an octet 
baryon. The lowest decay modes are: A — > N{p) + ir(— p), S* — > A(p) + n(—p), S* — > E(p) + ir(— p), as shown in 
Figure [8] Overall, most of the particles are stable. 

Finally, we measure the renormalized fermion anisotropy on the Nf = 2+1 lattices. We tuned the fermion anisotropy 

in the three-flavor calculation in Ref. [Til ] , where we found that the fermion action coefficients are consistent for bare 

PCAC quark masses up to about 175 MeV. Figure [5] shows the meson dispersion on the 24 3 x 128, atnii = —0.0840 and 

a t m s = —0.0743 ensembles. The effective mass plots are shown for the ground-state principal correlators at momenta 

p = j^-n with n £ {0, 1, 2, 3}; the fitted range and extracted energies are shown as straight lines across the effective 

mass plots. The inset shows the fitted renormalized fermion anisotropy at each n 2 . The speed of light c is measured 

from the energy of the boosted hadron using aiEHip) 2 = ci 2 Eh(0) 2 + jr^rn 2 . The values of c from the pion and 

. R 

rho mesons are about two and three standard deviations away from unity. Such a small deviation is also expected 
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FIG. 7: All measured meson masses as functions of the squared pseudoscalar masses. The diamonds and squares are measured 
with m s = —0.0540 but with two different volumes, 12 3 x 96 and 16 3 x 96; the upward-pointing triangles are those with 
m s = —0.0618 and 12 3 x 96 volume; the downward triangles and pentagons are measured with m s = —0.0743 and two different 
volumes, 16 3 x 128 and 24 3 x 128. The (red) dot-dashed lines indicate the decay thresholds for the 12 3 (upper) and 16 3 (lower) 
m s = —0.0540 ensembles, while the (blue) dashed lines are for the 16 3 (upper) and 24 3 (lower) m B = —0.0743. The lowest 
decay thresholds are: p — > 7r(p) + 7r(— p), K* — > K(p) + 7r(— p), (j> — » K(p) + K(—p), ao —> 7r(0) + rj(0), ai — > 7r(0) + p(0) and 
bi — > 7r(0) + uj(0) (with uj approximated by p) where the minimum allowed momentum p on the lattice is |p. 



on isotropic lattices. For example, c for the pion and rho are about two standard deviations away from unity on the 
MILC coarse asqtad lattice ensembles [29j. 



B. Static-Quark Potential 

V(r), the energy of two static color sources separated by distance r provides a useful reference scale for spectrum 
calculations. This is most usefully described by the Sommer parameter ro, defined by the condition 

-r 2 ^| r=ro = 1.65. (19) 

The potential is computed by measuring correlations between operators creating a static color source in the funda- 
mental representation of SU(3), connected via a gauge covariant parallel tranporter to a source in the 3 representation. 
The gauge connector can be formed by any sum of path-ordered products of link variables that respects the symme- 
try of rotations about the inter-source axis. Better ground-state operators are formed by using stout-smeared link 
variables in the path-ordered connections and by using an operator optimized using the variational method. 

A basis of five operators is constructed from the set of straight connectors and staples linking the mid-point 
between the two color sources. In forming the temporal correlators, straight, unsmcarcd temporal links arc used for 
the propagator of the static source. A five-by-five correlation matrix, Gij(r,t) is then computed for a range of time 
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FIG. 8: All measured baryon masses as functions of the squared pseudoscalar masses. The diamonds and squares are measured 
with m B = —0.0540 but with two different volumes, 12 3 x 96 and 16 3 x 96; the upward-pointing triangles are those with 
m s = —0.0618 and 12 3 x 96 volume; the downward triangles and pentagons are measured with m s = —0.0743 and two different 
volumes, 16 3 x 128 and 24 3 x 128. The (red) dot-dashed lines indicate the decay thresholds for the 12 3 (upper) and 16 3 (lower) 
m s = —0.0540 ensembles, while the (blue) dashed lines are for the 16 3 (upper) and 24 3 (lower) m B 



decay threshold are: A — > N(p) + n(—p), E* — ► A(p) + n(- 
on the lattice is j^-. 



-0.0743. The lowest 
H(p) + 7r(— p) where the minimum allowed momentum p 



separations t and values of r <E {1, N x /2} along a lattice axis. As outlined in the previous sub-section, this correlation 
matrix can be analysed using the variational method to make a more reliable ground-state energy extraction. 

Once the potential energy for a range of values of r has been determined, the data are compared with the Cornell 
model, 

V(r) = V + - + or, (20) 
r 

and best-fit values for the parameters a, a and Vq arc determined. See an example from one of our ensembles 
{atmi = —0.0808 and atm s = —0.0743) in Figure [TOl In all cases, a small range of r values that span rg are used. 
Once values of these parameters are computed, a value of r was derived from Eq. [151 The QCD flux-tube is expected 
to break in the presence of dynamical quarks and the ground-state of the system should not be modelled by the 
Cornell potential at large separations. We fit the data to Eq. [19] in a sufficiently small range of r such that this issue 
does not arise. No good evidence of this "string-breaking" effect was observed in our data at larger separations. This 
observation fits with previous investigations [30j, which established the need to include appropriate operators that 
construct two disconnected static-source-light-quark systems to measure the full spectrum. 
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FIG. 9: Effective-mass piots using the ground-state principai correiators from the pion and rho-meson masses at 4 different 
momenta on the iightest ensembie (m s = —0.0743) with voiume 24 3 x 128. The insets show the energy squared in temporai 
iattice units versus n 2 , which is reiated to the momentum by p 2 = ^-n 2 . 




FIG. 10: Resuits for the static-quark potentiai for the atirii = —0.0808, atm B — —0.0743 mass set. The ieft panel shows the 
effective energies, atV c s(r) for each r. The right panel shows the resulting fit to the potential using Eq. [20] 



IV. CHOOSING THE BARE STRANGE-QUARK MASS 

The appropriate value for the strange-quark mass in the lattice action is not known a priori. The Wilson formulation 
makes the task of choosing a sensible value for this parameter more difficult, as the breaking of chiral symmetry at 
the action level induces an additive mass renormalization. In dynamical simulations, changes to the strange-quark 
mass parameter in the action cause all observablcs to change. We suggest a helpful starting point for solving this issue 
is to determine where reference simulations lie in a parameterized two-dimensional coordinate system. Note that to 
leading order in chiral perturbation theory, the pseudoscalar masses are related to the quark masses via 

m% = 2B(m qi +m q2 ), (21) 

where B is a low-energy constant and m qi are the quark masses that compose the meson. The light-quark dependence 
can be eliminated using the linear combination (2m^- — m%). A useful property of a new coordinate system would 
be to remove all explicit dependence on the lattice cut-off. Such a dependence can be suppressed (if not completely 
removed) by taking ratios of hadron masses. One good candidate is the £1 baryon mass, which is stable against QCD 
decays and which has a simplified chiral extrapolation due to its lack of light valence quarks. An alternative is the 
S(l/2), which also decays only weakly and is statistically clean to measure. Appendix A shows a comparison between 
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-0.0826 
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3.221(25) 
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96 
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3.110(31) 
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96 
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3.647(14) 
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-0.0743 
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-0.0743 


-0.0743 


3.214(10) 



TABLE VIII: Sommer scale r /a B . 



N s 


N t 


a t mi 


a t m s 


In 


sn 




a t m,K 


a t mn 




12 


96 


-0.0540 


-0.0540 


0.564(14) 


0.564(14) 


0.2781(9) 


0.2781(9) 


0.556(7) 


0.334(3) 


12 


96 


-0.0699 


-0.0540 


0.356(8) 


0.535(10) 


0.1992(17) 


0.2227(15) 


0.501(4) 


0.3031(18) 


12 


96 


-0.0794 


-0.0540 


0.214(6) 


0.532(11) 


0.1393(17) 


0.1841(13) 


0.452(5) 


0.268(3) 


12 


96 


-0.0826 


-0.0540 


0.148(6) 


0.498(13) 


0.1144(19) 


0.1691(17) 


0.447(5) 
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16 
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0.161(9) 
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0.2322(15) 


0.2322(15) 
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0.286(5) 


16 
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-0.0743 


-0.0743 


0.397(7) 


0.397(7) 


0.1483(2) 


0.1483(2) 


0.353(3) 


0.2159(6) 
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-0.0808 


-0.0743 


0.231(6) 
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0.311(6) 


0.1901(18) 
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0.154(4) 
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0.1032(5) 


0.304(3) 


0.1845(11) 


16 


128 


-0.0840 


-0.0743 


0.124(4) 
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0.0691(6) 


0.0970(5) 


0.294(4) 


0.1788(13) 
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0.1205(15) 
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0.2945(16) 


0.1788(6) 



TABLE IX: Values of In and sn. 



these two choices. Therefore, we suggest two dimensionless coordinates, Iq and sq: 

9ml 



In = T-f, (22) 
9{2m\ - ml) 

SQ = ' (23) 

where the factor of 9/4 is a convenient normalization which makes Iq = sq = 1 in the static-quark limit. Note that 
three-flavor-degenerate theories lie on the diagonal line across the unit square. Table HXl summarizes all the {1,s}q 
values calculated in this work. Hadron masses are taken from Tables fVll and IVlTl in Sec. fVl 

In Figure [TT1 we locate all the simulations performed in this work using their Iq-sq coordinates. The dashed line 
runs horizontally from the physical point. We add two more strange-mass candidates: —0.0618, —0.0743 which are 
the points on the diagonal line. The choice of —0.0743 seems to anchor the correct m s value for Nf = 3 within one 
standard deviation of physical. Since we expect only a few percent deviation coming from the next-to-lcading effects 
on sq, we settle on atm s = —0.0743 for our final choice of strange bare mass; the points to the left of the Nf = 3 
points are Nf = 2 + 1 points with fixed strange input parameters. At the lightest simulation point, sq only differs 
from Nf = 3 by less than 2a. The running of the quantity sn is indeed small and is thus a good means for tuning 
the strange-quark mass for fixed-/3 simulation. We note however, that the trajectory followed by simulations as bare 
lattice parameters are changed is dependent on the details of the lattice action and is not universal; different actions 
may follow different paths as their bare parameters change. 
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FIG. 11: The location of the dynamical ensembles used in this work in the sn-ln plane. The circle (black) indicates the physical 
point s£ hys }. The (red) diamonds and squares are generated on 12 3 x 96 and 16 3 x 96 lattices with a t m s — —0.0540; 

the (green) upper triangle is the ensemble on a 12 3 x 96 lattice with atm B — —0.0618 and (blue) upside-down triangles and 



pentagons represent the ensembles on 16 x 128 and 24 x 128 lattices with atm s 



-0.0743. Detailed parameters can be found 



in in Table lIXl The horizontal dashed (pink) line indicates constant sn from the physical point and the diagonal line indicates 
three-flavor degenerate theories. 



V. EXTRAPOLATION TO THE PHYSICAL QUARK MASSES 



Following the discussion in Sec. IIV1 we adopt the coordinates Iq and so to perform extrapolation of the meson and 
baryon masses. To avoid the ambiguity in the lattice-spacing determination, we extrapolate mass ratios using 
the simplest ansatz consistent with leading order chiral effective theory, 



= c + ciln + c s sq 

tuq 



(24) 



with n = 2 for pseudoscalar mesons and n = 1 for all other hadrons. With such a parameterization, care is needed to 
take account of the statistical errors of Iq and so in the fit. Consider a general fit of the form / = a + bx + cy where 
/, x and y are all quantities with statistical error. We wish to find the combination of a, b, c which minimizes 



E 



(f(a,b,c;xi,yi) - (fr))' 



(25) 



where i indexes different data points {x,y, /}, (...) indicates a mean over all configurations and a is the statistical 
error of each quantity. The extrapolation (minimizing a quantity as in Eq. [25] with / = x = Iq and y = so) is 
taken to physical {I, s}o and we then take mo as experimental input to make physical predictions. 



Hadrons 



The % 2 /dof for the fits of hadronic data are all around or smaller than 1 for both meson and baryon masses. 
Figures [T2l and [T3l show the "sliced" plots of selected mass ratio with fixed Iq (or so)- The so are almost a constant 
for the same sea atm s ] this is why we see almost a single extrapolated line in the left column of the figures. The 
hadron masses linearly increase with Iq and decrease with so- The ratio of TOs/mo is almost constant with respect to 
so, indicating its insensitivity to the sea strange mass. The strange-mass dependence is almost completely canceled 
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all 



a t m s = -0.0743 



a t m,s = -0.0540 



m K 
m v 
m p 
ma* 
m<{, 
m ao 
m ai 



n/a 

0.570(5) [0.04] (4.29) 
0.780(8)[0.72](1.27) 
0.896(7) [0.49] (0.5) 
1.011(6)[0.42](0.84) 
0.98(6)[0.3](0.13) 
1.19(3)[0.7](3.39) 
1.26(3) [1.09] (2.28) 



0.476(6)[0.14](3.74) 
0.546(5)[0.25](0.2) 
0.812(8)[0.12](5.4) 
0.912(7) [0.1] (2.27) 
1.012(7)[0.15](0.75) 
0.98(7)[0.11](0.49) 

1.23(3)[1.1](0.1) 
1.33(4)[0.41](8.47) 



0.578(13)[0.17](16.7) 
0.677(12)[0.62](23.79) 
0.64(3) [2.96] (16.52) 
0.828(20)[1.7](7.19) 
1.007(18)[0.97](1.3) 
1.06(17)[0.05](7.73) 
1.03(7)[0.09](16.17) 
1.03(7)[0.53](16.18) 



ms 
m A 
mA 



1.020(12)[0.49](8.48) 
1.216(10)[0.62](2.15) 

1.319(9) [0.8] (0.08) 
1.166(10) [1.18] (4.51) 
1.325(12) [0.97] (7.57) 
1.461(9)[1.07](5.52) 

1.566(6)[1.](2.17) 



1.033(13)[0 
1.226(11)[0 
1.309(11)[0 
1.167(12)[1 
1.367(15)[0 
1.491(9)[0 
1.582(5)[0 



26](9.87) 
98] (3.03) 
63] (0.69) 
57](4.56) 
36](10.93) 
4] (7.67) 
84] (3.16) 



0.93(3) [0.43] (0.72) 
1.16(2) [0.5] (2.74) 

1.334(20)[0.62](1.22) 
1.13(2)[0.53](1.27) 
1.16(2) [4.81] (5.78) 
1.34(2)[2.47](3.06) 

1.506(18)[1.81](1.79) 



TABLE X: Hadron masses (in GeV) obtained from (m,H /mn) n (n = 2 for pseudoscalar mesons and 1 for the other hadrons) 
extrapolations using different sea-strange ensembles. The square brackets indicate the xVdof for each fit (with dof = 6), and 
the second parentheses indicate the percent deviation of the central value from experimental values. 



out in such a combination. Later in the appendix, we see the quantity (2m|- — m^.)/m| is relatively constant with 
respect to changes in m^. 

Figure [TJ] and the first column in Table [X] summarize all of our extrapolated masses along with the experimental 
values. The second half of the plots shows the relative discrepancy in percent between our calculation and the exper- 
imental numbers. The meson sector appears to have good agreement with experiment; overall, 0.1-4.3% discrepancy 
from experimental values. The biggest discrepancy comes in the rj, which we estimate using a combination of light 
and strange pesudoscalar mesons. All the vector mesons are in good consistency with experiment; the p meson is 
only 1.2cr away. These vector mesons are below the decay threshold on our ensemble; no decays are observed. The 
extrapolated scalar meson ao is consistent with the resonance at 980 MeV. The b\ meson is slightly higher than a±, 
and both of them are 2-3cr away from experiment. 

The baryon sector, on the other hand, does not work as well as the meson extrapolation. Non-strange baryons, 
such as nucleon and Delta, have the biggest discrepancy, by as much as 8.5er. This is likely due to contributions from 
next-to- leading-order chiral perturbation theory (or pion-loop contributions), which are not as negligible as the meson 
ones. This becomes evident as we increase the number of strange quarks in the baryon: the discrepancy is smaller 
in the Sigma and cascade. To have better control of the chiral extrapolation to higher order, we must have better 
statistics on these measurements; this will be a task for the near future once we complete all of our gauge generation. 

Finally, we compare the extrapolation results using all atm s ensembles and using a single ensemble of either 
atm s = —0.0540 or atm s = —0.0743 alone; results are summarized in Table Ixl In both cases, the 4> measurements 
are in good agreement with experiment; this is expected once f2 is fixed to 1.672 GeV. However, the kaon masses 
from atm s = —0.0540 ensemble are almost 17% away from experimental ones. This is also not surprising since the 
atm s = —0.0540 ensemble was selected using the J parameter strange-quark mass setting, where MILC had seen 
14-25% discrepancy in the strange-quark mass tuning. The kaon mass from the atm s = —0.0743 ensemble, on the 
other hand, is only 3<r away from the physical one, which is relatively close for a tuning using degenerate light and 
strange masses. The extrapolations using a t m s = —0.0743 alone versus all a t m s ensembles are in rough agreement 
within a few u, indicating that atm s = —0.0743 is a good candidate for gauge generation. 



B. Sommer Scale at the Physical Quark Masses 



Using the static-potential data in Table [VLLTl we can extrapolate r^ma to the physical limit using {/, s}n coordinates 
and the simplest functional form: 

romn = f + fiki + f s sn- (26) 
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FIG. 12: Selected meson mass ratios as functions of Iq and sq. Differently shaded (or colored) points correspond to the 
atmi tS combinations in Figure [Pol detailed numbers can be found in Table [VTl The smaller-volume ensembles {atmi, atm 3 } — 
{-0.0826, -0.0540} and {-0.0840, -0.0743} are excluded from the fits. The lines indicate the "projected" leading chiral 
extrapolation fit in Iq and sq while keeping the other one fixed. The black (circular) point is the extrapolated point at physical 
In and sq. 
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FIG. 13: Selected baryon mass ratios as functions of the Iq and sn- Differently shaded (or colored) points correspond to the 
atmi tS combinations in Figure [151 detailed numbers can be found in Table IVTT1 The smaller-volume ensembles {atmi, a t m s } = 
{-0.0826, -0.0540} and {-0.0840, -0.0743} are excluded from the fits. The lines indicate the "projected" leading chiral 
extrapolation fit in In and sq while keeping the other one fixed. The black (circular) point is the extrapolated point at physical 
In and sn- 
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FIG. 14: Summary of the extrapolated hadron masses compared with their experimental values. 
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FIG. 15: The assignment of colors from different ensembles in coordinates at{mi,m s }. The convention will remain consistent 
when used again for later extrapolations. 



Once this parameterization is known, too serves as experimental input and tq becomes a physical prediction. Using 
this technique, we determine the Sommer scale in the physical limit, rQ hys . We compose a dimensionless ratio 
^■(atrnfi) and extrapolate using Eq. 1251 After a chiral extrapolation including all three strange ensembles, we find 
the dimensionless ratio romji/^i? = 1.100(11). The "sliced" fits projected on to a single parameter Iq (left) and sq 
(right) at each point are shown in Figure [TBI Then we substitute in the physical Omega mass to find 

rf ys = 0.454(5) fm, (27) 

with % 2 /dof = 1.5(0.7) and dof = 6. The biggest x 2 contribution comes from the atm s = —0.0618 ensemble. If we 
drop it, we improve x 2 /dof to 0.92(0.60) but find only a small change to 7-Q hys = 0.451(5) fm. Extrapolating using a 
ratio with rather than too gives rQ hys = 0.446(4) fm. These data are consistent with the MILC result in Ref. [3lj | . 
which gave a continuum-extrapolated value of 0.462(12) fm. 



C. Scale Setting 

We determined the lattice cut-off (in physical units) at the physical point. Again, we follow the strategy of using 
the Omega-baryon to set the scale in physical units. A best fit of all simulation data to the model 

a t m n = d + diln + d s s n (28) 

has x 2 /Ndof = 3.10. Inputting the physical coordinates {?o,so} = {0.0153,0.379} yields a t = 0.03506(23) fm and 
a s = 0.1227(8) fm. The low quality of the fit provides us with further incentive to avoid expressing the lattice cut-off 
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FIG. 16: "Slices" of romn/^R as functions of In (left) and sn (middle). The straight lines are the fitted functions according to 
Eq. 1251 The assignment of colors is shown in the Figure 1151 

scale in physical units except in an extrapolation to the continuum. No continuum extrapolation is possible with our 
current data set, since all our ensembles have a common value of the gauge coupling (3. 

VI. CONCLUSION AND OUTLOOK 

This paper presents our first investigation of a number of states in the light hadron spectrum of QCD with Nj = 2+1 
dynamical flavors. Simulations were performed on anisotropic lattices with the ratio of spatial and temporal scales 
fixed non-pcrturbatively to a s /a t = 3.5. 

The focus of this work has been to test a simple method for determining the bare strange-quark mass, to allow us to 
approach the physical theory. Conventionally, this has been difficult to achieve, as there is delicate coupling between 
lattice action parameters and the cut-off scale. We found it extremely useful to introduce a pair of coordinates, sq and 
In to parameterize the two-dimensional space of quark mass values. The degenerate three-flavor theory corresponds 
to the line Iq = sq. To leading order in the chiral effective theory, these two coordinates are proportional to the 
strange and light-quark masses. These coordinates have been shown to be useful for our simulations since sq shows 
mild dependence on changes to the lattice light-quark mass. This shows that a good approximate value of the strange 
quark mass can be found by following the three-flavor degenerate line to the point where sn takes its physical value 
before changing the light-quark masses, a strategy adopted in this calculation. 

With a lattice strange-quark mass close to the physical value, a number of the simplest light hadrons were inves- 
tigated. The finite-volume effects for the mesons were checked on our data sets and found to be mild. Of the states 
we investigated, the ao,a\ and b\ mesons could have decayed on our lattices. However, checking the overlap factors 
of the interpolating operators with the ground states on different volumes suggests that the states we measured are 
predominantly resonances. For the octet and decuplet baryons, a similar analysis predicts that the A, E* and S* 
are stable in our study. At the heavy strange-quark masses where calculations on 12 3 lattices were performed, some 
finite-volume effects were seen, but they were negligible on the larger lattice volumes. 

Physical predictions have been made by extrapolating simulation data as a function of these coordinates to the 
physical point, {Iq, sq} = {0.0153, 0.379}. These extrapolations have been seen to be robust and have the advantage 
of making no reference to the lattice cut-off. This should enable reliable contact with chiral effective theories to 
be made. In this analysis, only the most naive extrapolations have been performed, and some discrepancy between 
extrapolated hadron masses and experimental data remains. It is encouraging to note that at worst, this discrepancy 
is less than 5% for mesons. The largest mismatch occurs in the nucleon-mass determination, which disagrees with 
experiment by 8%. It is very likely that the use of a naive extrapolation is responsible. No extrapolation to the 
continuum limit has been carried out; at present, calculations at a single value of the gauge coupling f3 have been 
performed, so no such analysis is possible. 

The collaboration has begun to explore more challenging measurements on the ensembles described in this work. 
The anisotropic lattice should allow us to resolve heavier excited states and those states which have traditionally been 
statistically rather imprecise with better accuracy. These more difficult calculations include the hybrid and isoscalar 
mesons, including the glucballs. We are confident that a detailed picture of a broad range of light-hadron physics will 
emerge soon from these analyses. 
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APPENDIX A: ALTERNATIVE POSSIBILITIES FOR {l,s} x 

In this work, we have been using the dimensionlcss parameters {Z,s}n (defined in Eg. I22[) to set the strange-quark 
mass and to extrapolate hadron mass ratios. In this section, we discuss alternatives to the Q: 1. the strange vector 
meson 4>, which like the fl contains no valence up or down quarks; 2. the octet S which is statistically cleaner to 
measure than the decuplct f2; 3. a linear combination of octet baryons, 2S — S, where the linear combination is 
selected to hopefully cancel out the leading-order up/down-quark dependence. 

We first look at the sx dependence on the sea strange mass, a similar strategy as described in Sec.|lVl Figure [T7l is 
a similar to Figure QTJ which we used to tune the sea strange quarks. The blue up-side-down triangles (V = 16 3 ) and 
pentagon (24 3 ) are points from the atm s = —0.0743 ensembles, the red diamonds (V = 12 3 ) and squares (16 3 ) are 
from atm s = —0.0540 and the green triangles are from atm s = —0.0618. The leftmost plot corresponds to X = 4>, and 
it shows the strong dependence of sx on a t m s that we are looking for to tune the strange-quark mass. The X = 4> 
could be an alternative for setting the strange quark mass. This result also agrees with our choice of sea strange at 
a t m s = —0.0743. The results from sh (middle plot) show negligible sea strange dependence on Nf = 2 + 1; this makes 
it a poor index for tuning the bare strange quark, since we cannot distinguish when sea quarks are not degenerate 
anymore. Somehow the remaining freedom of the light quark in the S baryon dominates the chiral behavior, since 
the atm s = —0.0743 set is running up toward the physical line. This running is greatly improved when we consider 
X = 2S — S, which lies on the physical line for all lx- However, it is a poor candidate for tuning, since it shows little 
dependence on sx- 

Let us move on to how the extrapolation behavior depends on the choice of X. Table IXll summarizes the results for 
all choices of X. The extrapolation (performed according to the minimization process described in Eq. I25p uses all 
three atm s ensembles without the smaller volume on the lightest ensemble of atm s — —0.0540 and atm s = —0.0743 
sets. The X = 2S — E has the poorest x 2 /dof among them all; we will throw it away for reliable extrapolation 
comparison. The X = 3 fits have similar % 2 /dof to the <p but slightly worse. This is possibly due to its insensitivity 
to sh during the extrapolation. The X = 4> should be quantitatively comparable to X = Q coordinates. However, due 
to the lightness of the 4> mass, it is not difficult to see that the next-leading-order contributions to the extrapolation 
form would be larger than the X — ft, causing it to be a slightly poorer fit at leading order. Even though the fit 
using X — 4> has smaller statistical error, we expect the systematic error to be higher than X = O. We will leave 
estimation of the systcmatics to future precision calculations where statistical error will be more reasonable. Still, we 
see good consistency between X = 4> and X = Q results, which reinforces our belief in the stability of extrapolations 
using the dimensionless coordinates {/, s}n- 
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FIG. 17: The s x -lx plot with X = (f> (left) 3 (middle) and the linear combination 2H - E for Nf — 3 and Nf = 2 + 1 at /3 = 1.5. 
The circle (black) indicates the physical point {^ ya , s^ ys }. The (red) diamonds and squares are 12 3 x 96 and 16 3 x 96 from 
atm 3 = —0.0540 ensembles; (green) upper triangles are from 12 3 x 96 from atm 3 = —0.0618 and (blue) upside-down triangles 
and pentagons are from 16 3 x 128 and 24 3 x 128 from a t m a = —0.0743 ensembles. Detailed parameters can be found in Table ITXl 
The horizontal dashed (pink) line indicates physical sx, and the straight diagonal lines indicate the SU(3) limit. 
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TABLE XI: Hadron masses (in GeV) obtained from (ran /mx)" (n = 2 for pseudoscalar mesons and 1 for the other hadrons) 
extrapolations in terms of {I, s}x with X £ {fi, 4>, 3, 2H — S} using all sea-strange ensembles. The square brackets indicate the 
X 2 /dof on the fit and the second parentheses denote the central value deviations from experimental values in percent. 



atm v 
atfUp 
atrriK* 
atm,/, 
atm ao 
a t m ai 
atirib 1 



a t m p 
atms 
a t ms 
atm a 
OirriA 

atms* 
a t ma 



x = n 



x 



X = s 



X = 23 



0.570(5) [0.04] (4.29) 
0.780(8) [0.72] (1.27) 
0.896(7) [0.49] (0.5) 
1.011(6)[0.42](0.84) 
0.98(6)[0.3](0.13) 
1.19(3) [0.7] (3.39) 
1.26(3)[1.39](2.28) 



0.569(2) [0.14] (4.04) 
0.795(5)[2.08](3.21) 
0.907(3) [1.77] (1.69) 
n/a 

0.99(6) [0.29] (0.35) 
1.21(3)[0.76](2.02) 
1.28(3)[1.85](4.06) 



0.570(3)[0.1](4.2) 
0.785(7) [2.44] (1.93) 
0.905(5)[2.3](1.42) 
1.023(5)[1.92](0.32) 
0.97(7)[0.33](1.98) 
1.18(3)[1.42](3.88) 

1.27(3) [1.9] (2.89) 



0.569(3)[0.11](4.1) 
0.819(7)[3.95](6.31) 
0.933(5)[3.73](4.55) 
1.044(5) [3.] (2.32) 
1.00(7)[0.31](1.6) 
1.23(3)[1.73](0.32) 
1.32(3)[2.58](7.03) 



1.020(12)[0.49](8.48) 
1.216(10)[0.62](2.15) 

1.319(9)[0.8](0.08) 
1.166(10) [1.18] (4.51) 
1.325(12)[0.97](7.57) 
1.461(9) [1.07] (5.52) 

1.566(6)[1.](2.17) 
n/a 



1.029(10)[0.89](9.45) 
1.226(8)[1.21](3.06) 
1.323(7) [1.42] (0.39) 
1.176(8)[1.86](5.38) 
1.335(17)[1.38](8.36) 
1.464(15) [1.41](5.72) 
1.568(12) [1.03] (2.3) 
1.685(10)[0.37](0.8) 



1.007(6) [0.06] (7.09) 
1.211(4) [2.7] (1.75) 
n/a 

1.161(3)[0.54](4.06) 
1.312(18)[1.68](6.51) 
1.450(15) [1.41] (4.66) 
1.561(13)[1.15](1.8) 
1.683(12) [0.89] (0.64) 



1.048(9) [0 
1.243(7) [3 
1.345(4) [3 
1.194(6)[1 
1.356(19)[2. 
1.490(17) [2 
1.593(15)[2 
1.708(12)[1 



63](11.5) 
81] (4.43) 
47] (2.02) 
6] (6.97) 
36] (10.07) 
49] (7.58) 
22] (3.94) 
82](2.17) 



APPENDIX B: STRANGE-SETTING COMPARISONS 

The J parameter is one common way to set the strange-quark mass; it is defined as 

dm v mK*(mj,-m p ) . . 

amp l{m K — m^) 

Here we examine how the parameter works for setting the strange-quark mass in our calculation. The upper four 
points in Figure [TBI are from N f = 2 + 1 at fixed m s = -0.0540 and mi G {-0.0699, -0.0794, -0.0826} from right to 
left with two volumes on —0.0826 (12 3 and 16 3 ). We hit the experimental value J cxp with the first mi = —0.0699, and 
the remaining points are within lc of J oxp . However, when we tried to extrapolate the kaon mass (see Sec. [VJ, we 
found that it missed the experimental value by about 17%. Such a mismatch resulting from tuning the strange mass 
using the J parameter has previously been reported in the literature. For example, MILC also found their lattice 
J parameter on their coarse and fine lattices agreed with J cxp , but after extrapolation they found the sea strange- 
quark mass to be off by 25% and 14% on the coarse and fine lattices respectively[3lj]. Although the discrepancy 
seems to become smaller for finer lattices, the J parameter does not seem to be an ideal quantity for strange-quark 
tuning. Similar conclusions can also be reached by observing the 4 lower points in Figure [T5] which correspond to a 
m s = —0.0743, Nf = 2 + 1 simulation; these are only two a away from the points for m s = —0.0540. The J parameter 
is not sensitive enough to changes in the strange sea-quark mass. 
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FIG. 18: J-parameter plot for Nj = 2 + 1 at /3 = 1.5. The upper diamond, upside-down triangle, pentagon (V = 12 3 x 96) 
and square (V = 16 3 x 96) points are from m s = —0.0540 ensembles and the triangle, lower diamond, upside-down triangle 
(V = 16 3 x 128) and square (V = 24 3 x 128) points are from m s = —0.0743; the circle (black) indicates the physical point 
{Zq ys , J phys }; the dashed line indicates the physical J value. 



